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NOTATION 

Lift  coefficient 
Free  stream  pressure 
Local  blade  radius 
Propeller  radius 
Time 

Velocity 

Circumferential  velocity 

Tangential  velocity  component  of  v/ake 

Axial  velocity  component  of  wake 

Free  stream  velocity 

Propeller  shaft  inclination  angle 

Blade  angular  position 

Cavitation  number 

Frequency  or  rotation 

Blade  rate  of  rotation 
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1 . 0  INTRODUCTION 

This  report  describes  work  carried  out  by  HYDRONAUTICS, 
Incorporated  as  part  of  a  study  of  partial  chord-length  cavities 
under  Office  of  Naval  Research  Contract  N00014-78-C-0144 .  This 
work  is  sponsored  under  the  U.S.  Navy  General  Hydromechanics 
Research  (GHR)  program  administered  by  the  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center  (DTNSRDC) . 

The  total  period  of  this  contract  has  been  four  years, 
beginning  in  April  1978  and  ending  in  Nay  1982.  During  the 
first  two  years  of  this  contract  a  theory  was  developed  at 
HYDRONAUTICS  giving  a  mathematical  representation  of  the  phe¬ 
nomenon  of  inception  of  leading-edge  cavity  as  well  as  its 
extent  on  hydrofoils  and  wings  Ln  uniform  flow.  During 
this  period  numerical  methods  were  developed,  based  on  this 
theory,  for  the  computation  of  leading-edge  cavity  on  an  air¬ 
foil  of  given  section  geometry  operating  in  uniform  flow  at 
zero  on  finite  angle  of  attack.  Computational  results  showed 
very  good  agreement  with  experiment.  References  1  and  5  give 
the  details  of  theory  development  and  the  correlation  of  the 
numerical  results  with  experimental  data. 

During  the  third  year  the  HYDRONAUTICS  cavity  flow  theory 
was  extended  to  permit  its  use  for  the  prediction  of  inception 
and  fluctuation  of  leading-edge  cavity  on  a  marine  propeller 
blade.  Using  the  theory  in  this  extended  form,  a  numerical 
prediction  was  made  of  the  inception  and  fluctuation  of  leading- 
edge  cavity  on  a  three-dimensional  marine  propeller  blade  under¬ 
going  a  purely  sinusoidal  fluctuation  of  angle  of  attack,  such 
as  occurs  with  operation  of  a  marine  propeller  with  an  inclined 
shaft,  and  uniform  inflow. 
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In  addition,  during  the  third  year,  the  cavity- flow  theory 
was  also  extended  to  make  it  applicable  to  noncavitating  flows 
involving  wakes  and  separation  bubbles.  A  paper  describing 
this  work  was  presented  in  October  1980  at  the  Thriteenth 
Symposium  on  Naval  Hydrodynamics  in  Tokyo,  Japan.  Details  of 
the  numerical  example  for  the  inclined-shaft  marine  propeller 
and  a  copy  of  this  Symposium  paper  were  included  in  Reference  6. 

During  the  final  year  of  this  contract,  a  validation  study 
for  a  numerical  method  for  predicting  unsteady,  partial  chord 
blade  cavities  on  marine  propellers,  which  is  described  in 
detail  in  this  report,  and  is  based  on  the  HYDRONAUTICS  cavity  flow 
theory,  was  conducted  for  typical  marine  propeller  operating  in  a 
realistic,  non-uniform  ship  wake.  Numerical  predictions  were 
made  for  the  inception  and  fluctuation  of  blade  leading  edge 
cavities  on  a  propeller  operating  in  the  wake  of  a  typical 
single -screw  cargo  ship.  The  predicted  results  were  then  com¬ 
pared  with  the  experimental  observation  of 
model  of  the  propeller  under  the  same  operating  condition, 

As  described  in  Section  4.0  of  this  report,  it  was  orig¬ 
inally  planned  to  consider  at  least  one  additional  propeller 
in  the  validation  study.  However,  there  was  not  found  another 
set  of  data  clearly  defining  suitable  leading  edge  cavitation 
over  a  significant  percent  of  blade  length  and  at  an  adequate 
number  of  angular  positions  of  the  blade.  An  inclined-shaft, 
uniform-flow  propeller  cavitation  test  was  then  conducted  in 
the  HYDRONAUTICS  high-speed,  free-surface  wave  channel  with  a 
container  ship  propelLer  model.  Unfortunately  suitable  leading 
edge  cavities  were  not  obtained  in  this  test. 

The  present  report  contains  the  results  of  numerical  pre¬ 
diction  of  leading-edge  cavity  made  for  the  inclined-shaft 
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propellcr  operating  in  uniform  flow  (as  presented  in  Reference  6) 
and  the  results  of  the  validation  study  made  for  the  cargo  ship 
prope 1 ler . 
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2.0  NUMERICAL  PREDICTION  OF  MARINE  PROPELLER  CAVITATION 

Numerical  predictions  were  made  of  the  inception  and  fluc¬ 
tuation  of  lending-edge  cavities  at  three  representative  radial 
blade  sections  of  the  selected  propellers.  Propeller  lifting- 
surface  theory  was  used,  in  conjunction  with  Lighthill's  Rule 
(Reference  1)  for  flow  at  foil  leading  edges,  to  predict  the 
distributions  of  perturbation  velocities  and  pressure  distribution 
due  to  blade  thickness  and  loading  for  each  of  the  selected 
blade  sections.  Calculations  were  carried  out  at  24  angular 
positions  of  the  propeller  blade,  using  even,  15 -degree  .incre¬ 
ments.  The  resulting,  velocity  distributions  were  then  used 
in  conjunction  with  the  blade  section  .1  ending -edge  cavity- flow 
theory  developed  at  HYDRONAUTICS  (References  1  and  5),  to  com¬ 
pute  unsteady  propeller  blade  cavity  lengths. 

2.1  Dose r L ption  of  Computational  Procedure 

For  each  selected  blade  section  and  each  of  the  twenty-four 
angular  positions  of  the  blade,  the  perturbation  velocities  due 
to  circulation  and  blade  thickness  were  obtained  using  computer 
program  PL'F-2,  Reference  3.  This  program  is  based  on  lifting- 
surface  theory  and  was  developed  at  MIT  for  the  computation  of 
unsteady  thrust  and  tongue  of  a  propeller  operating  in  non- 
uniform  flow.  For  each  selected  blade  section,  program  PUT- 2 
gave  the  perturbation  velocities  at  nine  chordwise  locations. 

A  continuous  chordwise  distribution  of  these  velocities  was 
obtained  from  the  nine  discreet  values  by  means  of  a  cubic 
spline  f i t . 

Of  the  nine  points  obtained  from  the  results  provided  by 
PUF-2,  the  point  nearest  to  the  leading  edge  was  2.5  percent  of 
chordlength  from  the  leading  edge.  A  distribution  of  the 
perturbation  velocities  was  needed  between  the  leading  edge 
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and  this  point.  This  distribution  was  obtained  using  Lighthill ' s 
Rule,  Reference  1.  In  order  to  use  I.ighthilL's  Rule  an  estimate 
of  the  blade  section  angle  of  attack  was  needed  at  each  angular 
position  of  the  blade.  This  estimate  was  made  in  the  manner 
described  below  for  each  blade  section  of  interest. 

Section  lift  coefficients  were  calculated  using,  a  computer 
program  based  on  the  lifting-] ino  theory  for  desired  value 
of  blade  section  and  blade  angular  position.  Using,  the  known 
physical  camber  of  each  blade  section  and  the  lifting-surface 
for  camber  presented  in  Reference  4,  the  effective  lift  coeffi¬ 
cient  due  to  camber  was  estimated.  The  section  lift  coefficient 
due  to  angle  of  attack  was  then  obtained  by  subtracting  this 
lift  coefficient  due  to  camber  from  the  total  section  lift  coef¬ 
ficient  determined  using  lifting-line  theory.  The  angle  of 
attack  for  the  equivalent  2-D  airfoil  section  was  then  obtained 
by  dividing  the  lift  coefficient  due  to  angle  of  attack  by  the 
section  lift- curve  slope. 

After  the  perturbation  velocities  were  obtained  Cor  each 
blade  section  and  blade  angular  positions,  computer  programs 
based  on  the  HYDRGNAUTICS  cavity-flow  theory,  as  described  in 
References  1  and  5,  were  used  to  predict  the  section  leading  edge 
cavity  geometry  at  each  section  and  angular  position.  For  each 
blade  section,  an  initial  estimate  of  cavity  volume  growth  (or 
collapse)  rate  was  made  by  the  quasi-steady  treatment  of  calcu¬ 
lating  the  cavity  volumes  at  the  twenty- four  angular  positions 
and  obtaining  the  local  rates  at  these  positions  by  using  the 
calculated  lengths.  These  growth  (or  collapse)  rates  of 
cavity-volume  were  then  incorporated  into  the  computations  and 
the  cavity  volumes  were  recalculated  at  the  twenty-four  angular 
positions.  The  new  rates  were  then  calculated.  This  procedure 
was  repeated  until  the  local  rate  of  growth  (or  collapse)  at 
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each  blade  position  converged  within  a  specified  tolerance. 

Corrected  cavity  growth  (or  collapse)  rates  due  to  fluctuation 
of  angle  of  attack  were  thus  obtained  at  the  selected  blade 
sections  . 

The  unsteady  cavity  flow  predictional  methods  which  have 
been  developed  as  part  of  this  contract ,  arc  two-dimensional 
methods.  It  has  been  previously  demonstrated  that  the  methods 
can  be  successfully  applied  to  finite  aspect  ratio  foils.  It 
was  hypothesized,  based  on  this  experience  with  three-dimensiona 
hydrotoils,  that  the  methods  should  be  applicable  to  marine 
propeller  blades  in  which  the  cavity  spanwisc  extent  was  signi¬ 
ficantly  greater  than  cavity  chordwise  extent  or  length  and  in 
which  this  cavity  length  did  not  change  too  rapidly  with  radius. 

The  method  is  only  applicable  to  leading  edge  cavities  which 
are  distinct  from  tip  vortex  cavities.  Validation  of  the 
applicability  of  these  cavity  flow  methods  to  t.he  prediction 
of  unsteady  propeller  cavitation  therefore  requires  test  data 
for  a  propeller  with  Leading  edge  cavities  of  this  type. 
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3.0  PREDICTION  OF  CAVITATION  ON 
MARINE  PROPELLER  ON  INCLINED  SHAFT 


Numerical  predictions  were  made  of  inception  and  time  varia¬ 
tion  of  leading-edge  cavitation  on  a  propeller  in  which  the 
blades  were  undergoing  a  purely  sinusoidal  fluctuation  of  angle 
of  attack.  This  was  done  by  considering  a  marine  propeller 
operating  on  an  inclined  shaft  in  uniform  flow.  The  computations 
were  carried  out  for  two  representative  blade  sections,  0,684R 
and  0.937R,  where  R  is  the  overall  radius  of  the  propeller. 

Using  the  methods  described  in  Section  2,  the  fluctuating  cavity 
lengths  and  volumes  were  obtained  at  these  two  blade  sections. 


3  .  1  Select  ion  o  f  Prope  Her 

A  typical  single-screw  cargo  ship  propeller  was  chosen  for 
this  study.  This  was  one  of  the  two  propellers  designed  by 
DTNSRDC,  for  which  model  tests  had  been  conducted  at  MIT  in  a 
realistic  axial-flow  ship  wake  in  order  to  study  their  cavita¬ 
tion  characteristics.  This  propeller  was  chosen  because  it 
was  also  planned  to  use  the  MIT  test  data  to  validate  the 
HYDRONAUTICS  cavity-flow  theory  for  the  case  with  a  typical 
axial  wake.  This  latter  effort  could  not  be  carried  out,  be¬ 
cause  these  data  could  not  be  obtained  from  MIT  (See  Section 
4.1).  The  principal  particulars  of  this  propeller  model  are 
given  below: 


Overall  Radius,  R 
Design  RPM,  N 
Rotational  Speed,  N 
No .  of  Blades 

Pitch-Diameter  Ratio  at  0.7  R 
Expanded  Area  Ratio 
Blade  Section  Type 
Blade  Camber  Type 


6  inches 
1500 

1500  rpm 
5 

1.077 
0.725 
NACA  66 

a  =  0 . 8  meanline 
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3 • 3  Description  of  Inflow 

The  orientation  of  propeller  with  respect  to  inflow  is 
illustrated  in  Figure  1.  The  flow  conditions  used  in  computations 
are  given  below: 


Frees tream  Flow  Velocity.  U 

*  CO 

Inclination  of  Flow  to  Propeller  Axis,  6 
Freestrenm  pressure,  p^ 

Kean  Cavitation  Number  at  0.684R 
Mccin  Cavitation  Number  at  0.93/R 


12 . 33  knots 
10  degrees 
11.8  psi 
0.67 
0.39 


As  shown  in  Figure  1,  the  axial  component  of  the  free- 
strcam  velocity,  V^,  is  constant  over  the  propeller  disk  and  is: 

V.  U  cos  6  . 

A  <» 

The  shaft  inclination  produces  a  tangential  or  circumferential 
velocity  component  which  has  a  mean  value  of  Dr,  where  Q  is 
the  constant  rotational  speed  of  the  propeller  and  r  the  radial 
distance  of  the  blade  section  from  the  axis  of  rotation.  This 
tangential  velocity  also  has  a  time-varying  component, 

Vrj.  sin  wt , 


which  is  independent  of  the  radial  location  of  the  blade  section. 
This  fluctuating  part  of  the  tangential  velocity  is  due  to  the 
propeller  axis  inclination  to  the  flow.  The  amplitude  and  phase 
ofthis  fluctating  velocity  component  are 

VT  =  sin  6 

and 

03 1  =  0,  respectively, 

where  0  is  the  blade  angular  position  as  defined  in  Figure  1, 

The  magnitude  of  the  resultant  inflow  velocity  at  a  blade 
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section  at  radius  r  is  then 


VR  - 


V 


V 


+ 


V  2 
c 


where 


V.  =  U  cos  >S 
A 

and 


V  =  Qr  +  V™  sin  ut 
c  T 

=  52r  -r  U  sini  sinO 

CO 


3.3  Calculated  Results 

The  predicted  results  are  given  in  Figures  2  through  6. 
Figure  2  shows  the  circumferential  fluctuation  of  lift  coeffi¬ 
cient  at  radii  0 . 684R  and  0.937R.  This  variation  is  due  to  the 
propeller  axis  inclination  to  the  uniform  flow,  resulting  in 
sinusoidal  variation  of  blade  section  geometric  angle  of  attack. 
The  resulting  cavity  length  variations  at  the  two  blade  sections 
arc  shown  in  Figures  3  and  4.  Figure  5  gives  the  angular  rate 
of  variation  of  cavity  volume  and  Figure  6  gives  the  fluctuation 
of  cavity  volume  with  respect  to  the  mean  volume. 

3 . 4  Discussion  of  Results 

As  seen  in  Figures  3  and  4,  the  propeller  blade  had  cavita¬ 
tion  at  the  two  radii  of  interest  at  all  angular  positions  of 
the  blade.  The  inception  (onset)  and  washoff  of  cavitation  due 
to  varying  angle  of  attack  therefore  were  not  seen  at  any 
blade  position,  because  the  blade  section  angle  of  attack  was 
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large  enough  at  all  positions  of  the  blade  to  produce  leading 
edge  cavitation.  As  seen  in  Figure  6,  the  fluctuation  of  cavity 
volume  or  cross-sectional  area  is  small  and  does  not  exceed  27 
percent  of  the  mean  volume  at  0.684R  or  30  percent  of  the  mean 
volume  at  0.937R. 
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4.0  VALIDATION  STUDY  FOR  PROPELLER 
IN  NONUNIFORM  FLOW 

L\ .  1  Selection  of  Propeller  and  Inflow 

The  selection  of  a  suitable  propeller  and  inflow  for  use 
in  validating  the  predictional  methods  developed  under  this  con¬ 
tract  was  largely  based  on  the  availability  of  a  set  of  experi¬ 
mental  data  which  gave  a  clear  definition  of  cavitation  on  pro¬ 
peller  blade  at  various  angular  positions  of  the  blade  while 
operating  in  a  well  defined  wake. 

At  the  time  this  validation  study  was  proposed,  it  was 
intended  that  the  study  be  conducted  for  two  configurations  of 
propeller  operation:  the  propeller  operating  on  an  inclined 
shaft,  in  a  horizontal  uniform  flow,  and  hence  undergoing  a  purely 
sinusoidal  variation  of  angle  of  attack  and  the  propeller  oper¬ 
ating  in  a  typical  single-screw  cargo  ship  wake,  which  was  dom¬ 
inated  by  variations  in  axial  velocity. 

A  search  of  the  available  literature  was  made  for  cavitation 
data  for  the  uni  form- flow,  inclined-shaft  propeller  operation; 
but  no  set  of  experimental  data  clearly  defining  the  cavitation 
on  propeller  blade  at  a  number  of  angular  positions  of  the  blade 
was  found.  An  inclined-shaft,  uniform- flow  propeller  eavitation 
test  was  then  conducted  in  the  HYDRONAUTICS  high-speed,  free- 
surface  water  channel  with  a  container  ship  propeller*  model  of 
HYDRONAUTICS  design.  Unfortunately,  no  suitable  data  were  ob¬ 
tained  from  these  model  test.  The  model  test  was  conducted  for 
a  shaft  inclination  of  7.5  degrees  and  for  various  combinations 
of  free-stream  velocity,  rpm  and  static  pressure.  Clear,  par¬ 
tial  leading-edge  cavitation  could  not  be  obtained  at  any  test 
condition.  Cavitation,  when  it  occurred,  remained  mainly  in 
the  tip  region,  distorted  by  tip-vortex  cavitation,  and  the 
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cavity  in  that  region  extended  over  the  whole  chordlcngth. 


For  validation  Cor  a  propeller  operating  in  a  typical  single 
screw  ship  wake,  it  was  originally  proposed  to  use  data  from 
recent  cavitation  tests  conducted  at  MIT,  for  the  Maritime  Admin¬ 
istration,  with  a  typical  single-screw  cargo  ship  propeller. 
Unfortunately,  MIT  was  unable  to  make  available  suitable  test 
data  for  this  study.  After  a  search  of  the  published  test  data, 
those  presented  in  Reference  2  were  selected  for  this  study. 

These  data  are  from  a  propeller  cavitation  test  performed  at 
the  Norwegian  Ship  Model  Experiment  Tank  in  the  wake  of  a  typical 
single-screw  cargo  ship.  Details  of  the  propeller  geometry 
were  supplied  by  Dr.  Erling  lluse,  author  of  Reference  2. 

4.2  Description  of  S ej. c_c ted  _Pr o peller 


The  propeller  model  desig 
Ship  Model  Experiment  Tank  was 
principal  particulars  of  this 

Overall  Radius,  R 
rpm ,  N 

Rotational.  Speed, 

No.  of  Blades 
Pitch -Diameter  Ratio  at 
Expanded  Area  Ratio 
Blade  Section  Type 
Blade  Camber  Type 


rutted  as  P-548  by  the  Norwegian 
selected  for  this  study.  The 
propeller  model  arc  given  below: 

:  3.73  inches 

:  1140 

1 140  rpm 
:  4 

0.7R:  1.077 

:  0.738 

:  B-Series  (Wageningen) 

:  a  =  0.8  meanline. 


Some  geometric  properties  of  the  propeller  blade  geometry 
arc  given  in  Figure  7.  Details  of  propeller  geometry  are  not 
presented  at  the  request  of  Dr.  Huse. 


4.3  Description  of  Ship  Wake  and  Its  Effect  on  Cavitation 

A  velocity  contour  map  of  ship  wake  including  the  location 
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of  the  propeller  is  given  in  Figure  8.  This  wake  is  represented 
by  contours  of  constant  axial  inflow  velocity  made  dimensionless 
with  respect  to  ,  the  freestream  velocity.  The  inflow  to  the 
propeller  in  this  wake  is  assumed  to  be  purely  axial,  i.e.,  the 
flow  velocity  vector  is  parallel  to  the  axis  of  rotation  of  the 
propeller . 

The  physical  properties  of  the  propeller  inflow  in  the 
model  experiment  set-up  were  as  follows: 


Freestream  Velocity,  U  :  14.11  ft/s 

Freestream  Pressure,  pu  :  13.35  psi 

Mean  Cavitation  Number  at  0.684R:  0.53 

Mean  Cavitation  Number  at  0.8 11R:  0.39 

Mean  Cavitation  Number  at  0.937R:  0.30 


4 . 4  Calculated  and  _Kxper_imontal_  R£i>uJ_tj3 

Results  from  numerical  prediction  and  experiment;  are  given 
in  Figures  9  through  14.  Figure  9  shows  the  circumferential 
variation  of  blade  section  Lift  coefficient  at  the  three  radii 
of  interest.  Figure  10  gives  the  fluctuation  of  cavity  volume, 
with  respect  to  the  maximum  volume,  at  each  of  these  radii. 
Figures  11,  12,  and  13  show  the  circumferential  variation  of 
cavity-length,  as  obtained  from  both  the  numerical  prediction 
and  the  model  test,  at  0.684R,  0.811R  and  0.937R  respectively. 

As  described  in  Section  2.1,  the  three-dimensional  propeller 
blade  problem  was  reduced  to  one  in  two  dimensions,  where  the 
calculations  of  the  equivalent  2-D  section  angle  of  attack 
required  the  two-dimensional  lift-curve  slope.  For  the  numerical 
results  presented  in  Figures  11-13,  the  2-D  lift-curve  slope 
used  was  that  given  by  inviscid  theory,  viz.,  2?r .  However, 
experimental  data  indicate  that  it  is  more  appropriate  to  use 
a  value  of  about  5.73  (0.1  per  degree)  for  the  lift  curve  slope 
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of  typical  propeller  blade  sections.  A  comparison  was  made,  for 
the  blade  section  at  0.684R,  between  predictions  of  cavity  made 
using  these  two  lift-curve  slopes.  Figure  14  shows  this  compari¬ 
son.  As  seen  in  this  figure,  the  predicted  leading-edge  cavity 
changed  only  slightly  with  this  change  in  the  value  of  lift- 
curve  slope,  and  the  difference  is  not  considered  significant. 

Figure  15  shows  the  variation  of  cavity  length  with  axial 
velocity,  for  the  first  quadrant  of  propeller  blade  rotation. 

4.5  Discussion  of  Results 

The  purpose  of  that  phase  of  the  contract  discussed  here 
was  validation  of  the  applicability  of  HYDRONAUTICS  cavity-flow 
theory,  to  the  prediction  of  unsteady  marine  propeller  blade 
leading  ed.ge  cavitation.  The  procedure  used  for  propeller 
cavitation  prediction,  ns  described  in  Section  2.1,  involved 
five  step:;,  each  of  which  significantly  influences  the  prediction 
of  cavitation.  These  steps  are: 

(a)  computation  of  blade  surface  perturbation  velocities 
using  a  slightly  modified  version  of  lifting  surface  computer 
program  PDF- 2; 

(b)  computation  of  blade  section  lift  coefficient  at  each 
angular  position  of  blade  using  a  lifting-line  computer  pro¬ 
gram; 

(c)  use  of  tabulated  lifting-surface  data  to  calculate  the 
lift  coefficient  due  to  camber  for  an  equivalent  2-D  airfoil 
section ; 

(d)  use  of  a  2-D  lift-curve  slope  to  obtain  the  section 
angle  of  attack; 

(e)  use  of  the  HYDRONAUTICS  cavity-flow  theory  for  com¬ 
putation  of  unsteady  cavity  length  and  volume. 
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As  seen  in  Figures  11,  12,  and  13,  the  theory  agrees  quite 
well  with  experiment  as  far  as  the  extent  of  cavitation.  There 
exist,  however ,  some  discrepancies  in  the  point  of  inception  of 
ca . a  Cation  and  phase  of  maximum  cavitation: 

(a)  Near  the  30 0°  position  of  the  blade,  the  theory  pre¬ 
dicts  inception  of  cavitation  significantly  earlier  than  is 
observed  in  experiment . 

(b)  In  the  vicinity  of  the  180°  position  of  the  blade, 
the  theory  predicts  cavitation  while  the  experiment  shows  no 
evidence  of  cavitation. 

(c)  Tiie  experimental  observations  lag,  in  phase,  the 
predicted  occurrences  of  maximum  cavitation.  This  phase  lag 
is  especially  noticeable  at  0.684R  and  0.8 HR  radii,  where 
the  experiments  show'  maximum  cavitation  at  15°  and  about  40° 
blade -posit ions ,  respectively,  while  the  theory  predicts 
maximum  cavitation  at  0°  for  alL  radii. 

There  are  three  possible  major  sources  of  the  disagree¬ 
ments  described  above.  These  are:  factors  influencing  the 
experimental  results;  computational  methods  used  in  program 
PUF-2  to  account  for  the  unsteady  effects  on  the  blade  loads 
and  surface  perturbation  velocities;  and  the  cavity-flow 
theory  and  associated  computational  mehtod  used  to  determine 
the  time-history  of  cavity  growth  (or  collapse)  at  a  given 
blade  position. 

Near  the  300°  location  of  the  blade,  the  experiment  shows 
a  rather  significant  delay  in  the  observed  inception  of  cavita¬ 
tion  relative  to  theoretical  prediction.  However,  except  at 
the  0.937R  radius,  where  the  cavitation  is  expected  to  be 
significantly  influenced  by  three-dimensional  effects,  the 
theoretical  and  experimental  blade-position  for  cavity  wash- 
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off  are  in  good  agreement. 

In  cavitation-tunnel  tests,  a  delay  in  cavitation  Inception 
will  occur  if  the  air  content  in  the  tunnel  water  is  not  at  least 
saturated,  with  some  free  air  bubbles.  The  experiment  of  tin's 
study  was  conducted  in  a  closed  jet  water  tunnel,  and  in  such 
tunnels  the  air  content  in  water  i.s  typically  less  than  the 
saturated  value.  Therefore,  it  is  possible  that  the  delay  in 
cavitation  inception  shown  in  experiment,  compared  with  the 
theoretical  prediction,  is  due  at  least  in  part  to  a  low  air 
content  in  Lhe  tunnel.  A  low  air  content  might  also  explain 
the  absence  of  cavitation  in  experiment  in  the  vicinity  of  the 
ISO'3  position  of  the  Made,  where  the  theory  predicts  the 
occurrence  of  cavitation. 

The  second  possible  source  to  the  observed  discrepancies 
is  the  method  employed  in  program  l'UF-2  to  account  for  unsteady 
effects.  Figure  15  shows,  for  the  first  quadrant  of  propeller 
blade  rotation,  the  variation  of  predicted  cavity  length  with 
axial  velocity  in  ship’s  wake.  Figure  15  shows  that,  as  the 
blade  enters  regions  of  increasing  axial  velocity,  the  pre¬ 
dicted  cavity  length  decreases  monotonical ly  from  its  maximum 
at  the  0-dcgree  Location  of  blade.  The  computational  procedure 
predicts  the  occurrence  of  maximum  cavity  length  at  the  0-dcgree 
location,  where  the  axial  velocity  is  a  minimum.  This  prediction 
is  consistent  with  the  physical  phenomenon  in  the  absence  of 
unsteady  effects,  in  which  situation  the  maximum  cavity  length 
would  occur  at  the  0-degree  location,  the  location  of  the  maximum 
blade  section  inflow  angle  of  attack.  On  the  other  hand,  the 
presence  of  any  significant  unsteady  effects  would  cause  the 
blade  location  for  maximum  cavity  length  to  shift  from  the 
location  of  minimum  axial  velocity.  This  phase  difference 
between  minimum  axial  velocity  and  maximum  cavity  length  does 
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exist  in  the  experimental  results  (Figures  11,  12,  and  13),  sug¬ 
gesting,  therefore,  the  presence  of  the  effects  of  unsteadiness. 
Because  the  prediction  shows  no  such  difference  in  phase,  it  is 
possible  that  the  unsteady  effects  were  not  adequately  taken 
into  account  by  program  PUF-2  in  the  computation  of  blade  sur¬ 
face  per truba t ion  velocities. 

It.  appears  that  no  validation  of  program  PUF-2  based  on 
experimentally  measured  forces  on  a  single  propeller  blade  as  a 
function  of  angular  positions  in  an  axial  wake  has  been  carried 
out.  A  validation  study  was  conducted  for  PUF-2  by  Boswell, 
et  al .  (Reference  7),  but  this  was  for  total  forces  on  the 
propeller  and  therefore  the  conclusions  derived  in  Reference  7 
do  not  apply  to  this  study. 

Finally,  the  method  used  in  computing  propeller  blade 
leading-edge  cavity  characteristics  may  have  contributed  to  the 
disagreement,  between  theory  and  experiment,  in  the  inception 
and  phase.  As  described  in  Section  2.1,  the  time-history  of 
cavity  growth  (or  collapse)  was  taken  into  account  in  the 
computation  through  an  iterative  procedure  for  the  convergence 
of  time -rate  of  cavity  fluctuation  at  a  given  blade  position. 
There  is  a  possibility  that  the  time  -  dependence  of  the  fluc¬ 
tuation  of  cavity  was  not  fully  accounted  for  by  this  procedure. 

The  possible  major  sources  of  the  disagreements  have  been 
discussed  above.  There  is  no  means  with  available  data  of 
determining  which  of  these  sources  is  (are)  the  primary  con¬ 
tributor  (s)  to  these  disagreements.  Therefore  a  more  definitive 
explanation  for  these  discrepancies  based  on  these  results  is 
not  presently  possible. 
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4 . 6  Cone  lus  ions  and  Roconimenda  t  ions 

Based  on  the  validation  effort  described  in  this  section, 
it  can  be  concluded  that: 

•  The  proposed  method  of  analysis  for  unsteady  propeller 
blade  leading  edge  cavitation  appears  capable  of 
predicting  such  cavity  flows  with  generally  good  accuracy. 

9  The  primary  differences  between  the  measured  and  pre¬ 
dicted  cavity  flow  fluctuations  are  in  the  positions 
of  cavity  onset  and  washoff,  and  in  the  absence  of  a 
small  region  of  observed  cavitation  when  the  blade  is 
in  the  vertical,  downward  position  (■  --  180'1). 

®  There  are  a  number  of  possible  causes  for  the  observed 
differences  between  predicted  and  measured  results, 
including  possible  scale  effects  in  the  model  test 
results  used. 

«  The  method  appears  to  he  a  uscfuL  engineering  tool 
for  predicting  the  inception,  extent  and  volume  of 
unsteady  blade,  leading  edge  cavities. 

In  order  to  further  improve  confidence  in  the  method, 
additional,  validation  is  needed.  In  particular,  validation 
based  on  tests  with  the  propeller  inclined  to  a  uniform 
inflow  are  desirable,  as  these  eliminate  the  effect  of  factors 
such  as  flow  shear  and  turbulence  which  are  present  with 
typical  nonuniform  wakes.  A  comparison  with  data  from  the 
carefully  conducted  experiments  at  MIT  would  be  very  useful. 
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ORIENTATION  OF  PROPELLER 


WITH  RCSPECT  TO  INFLOW 


VA  =  U  COS  5 
A  « 
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T  •*> 

RESOLUTION  OF  INFLOW  VELOCITY  VECTOR 

0° 


DIRECTION  OF  PROPELLER  ROTATION 

(AS  VIEWED  LOOKING  IN  THE  DIRECTION 
OPPOSINC  AXIAL  INFLOW,  VA) 


FIGURE  1  -  DESCRIPTION  OF  INFLOW  AND  PROPELLER  OPERATION 

(INCLINED-SHAFT  PROPELLER  CAVITATION  PREDICTION) 


I 


HYDRONAUTICS,  INCORPORATED 


AT  r/R  =  0.684: 

CL  -  0.2207 
mean 

EXTREME  FLUCTUATIONS  ARE 

9. 75%  (+) 

AND  9.50%  (-) 

CORRESPONDING  FLUCTUATIONS 
OF  ANGLE  OF  ATTACK  ARE 

11.25%  (+) 

AND  10.90%  (-) 


AT  r/R  =  0.937: 

C  =  0.1689 

mean 

EXTREME  FLUCTUATIONS  ARE 

9.25%  (+) 

AND  9.00%  (-) 

CORRESPONDING  FLUCTUATIONS 
OF  ANGLE  OF  ATTACK  ARE 

12.50%  (+) 

AND  12.10%  (-} 


FIGURE  2  -  VARIATION  OF  LIFT  COEFFICIENT 
(INCLINED-SHAFT  PROPELLER 
CAVITATION  PREDICTION) 


CAVITY  LENGTH /CHORD  LENGTH 
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FIGURE  3  -  FLUCTUATION  OF  CAVITY  LENGTH  AT  0.684R 
(INCLINED-SHAFT  PROPELLER  CAVITATION 
PREDICTION) 
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FIGURE  4  -  FLUCTUATION  OF  CAVITY  LENGTH  AT  0.937R 
(INCLINED-SHAFT  PROPELLER  CAVITATION 
PREDICTION) 
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FIGURE  5  -  RATE  OF  CHANGE  OF  CAVITY  VOLUME 

(INCLINED-SHAFT  PROPELLER  CAVITATION 
PREDICTION) 
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a  =  CAVITATION  NUMBER 
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FIGURE  G  -  FLUCTUATION  OF  CAVITY  VOLUME 

(INCLINED  SHAFT  PROPELLER  CAVITATION 
PREDICTION) 


FIGURE  7  -  SOME  GEOMETRIC  DETAILS  OF  PROPELLER  OPERATING  IN  SHIP  WAKE 
(TAKEN  FROM  REFERENCE  2} 


LIFT  COEFFICIENT, 
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FIGURE  9  -  VARIATION  OF  BLADE  SECTION  LIFT  COEFFICIENT 
(PROPELLER  IN  SNIP  WAKE) 


CAVITY  VOLUME  FLUCTUATION 
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FIGURE  10  -  PREDICTED  FLUCTUATION  OF  CAVITY  VOLUME 
(PROPELLER  IN  SHIP  WAKE) 
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FIGURE  11  -  FLUCTUATION  OF  CAVITY  LENGTH  AT  0.684  R 
(PROPELLER  IN  SHIP  WAKE) 
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FIGURE  12  -  FLUCTUATION  OF  CAVITY  LENGTH  AT  0.811R 
(PROPELLER  IN  SHIP  WAKE) 
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FIGURE  13  -  FLUCTUATION  OF  CAVITY  LENGTH  AT  0.937  R 
(PROPELLER  IN  SHIP  WAKE) 
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